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Abstract 

We consider a monostable time-delayed reaction-diffusion equation 
arising from population dynamics models. We let a small parameter tend 
to zero and investigate the behavior of the solutions. We construct ac- 
curate lower barriers — by using a non standard bistable approximation 
of the monostable problem — and upper barriers. As a consequence, we 
prove the convergence to a propagating interface. 

Key Words: time-delayed reaction-diffusion equation, delay differential 
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1 Introduction 

In this work we investigate the singular limit, as e — > + , of u e : [— et, oo) : 
R the solution of the delayed reaction-diffusion equation 



d t u(t,x) = eAu(t,x) + -[f(u{t-eT,x))-u(t,x)}, t > 0, x £ R N , (1.1) 
supplemented with the initial data of delayed type 

u (9, x) = (p(^,xj , -ST < 9 < 0, x £ R N . (1.2) 

Here r > is a given delay parameter; / : [0, oo) — > [0, oo) is a given increasing 
and monostable nonlinearity — see (jl.4p for precise assumptions; the initial 
data (f : [— r, 0] x R w — > R is a given smooth function — see Assumption ll.il 

Equation (|1.1[) is widely used in population dynamics models. In this con- 
text, u(t,x) denotes the density of individuals at time t and spatial location x. 



1 AMS Subject Classifications: 35K57, 35R10, 92D25. 

The authors are supported by the French Agence Nationale de la Recherche within the 
project IDEE (ANR-2010-0112-01). 
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The function / is the birth rate of the population. Note that the birth feed- 
back appears with some time delay in order to take into account the period of 
maturation to become adult. Finally the term — u corresponds to a normalized 
death rate, while e > is a scaling parameter. 

When / takes the form of the so-called Ricker's function 

f(u) = aue~ u , a > 1, (1.3) 

equation (|1.1[) is commonly referred as the diffusive Nicholson's blowflies equa- 
tion. This kind of equation has been intensively studied in the literature. The 
purely reactive part, namely the underlying delay differential equation, has at- 
tracted the attention of many researchers during the past decades (see for in- 
stance [TS] and references cited therein). On the other hand, the diffusive equa- 
tion has also been extensively investigated from the spatial propagation point 
of view, that is speed of spread, travelling wave solutions (we refer for instance 
to So and Zou [2D], So, Wu and Zou [ID], Thieme and Zhao [22], Fang and Zhao 
[7], and the references therein). 

In this work, we consider the monostable equation (jl.ip in the so-called 
monotonic regime. Precisely we assume that / : [0, oo) — > [0, oo) is a function 
of the class C 2 such that 

7(o) = o, /(i) = i, /'(o)>i, /'(i)<i, 

< /'(u)>0, V«G (0,1), (1.4) 
J(u)>u, Vug (0,1). 

In particular, u = and u = 1 solve (|1.1[> . If we come back to example (|1.3|) . 
assuming a G (1, e) implies that / satisfies (jl.4[) . with In a playing the role of 1. 

Let us observe that, when r = 0, equation (| 1 . 1 1) reduces to the monostable 
reaction-diffusion equation 

d t u(t, x) = eAu{t, x) + (u(t, x)) , (1.5) 

with F(u) := f(u)—u. In view of (jl.4p . the nonlinearity F exhibits a monostable 
dynamics, namely F(0) = F(l) = 0, F(u) > for all u G (0, 1), and F'(0) > 
while F'(l) < 0. Under these assumptions, solutions of (|1.5[) with compactly 
supported initial data have been considered first by Freidlin [9] with probabilis- 
tic tools, then by Evans and Souganidis [6] with Hamilton- Jacobi techniques 
(we also refer to [U [5] and the references therein). This problem has been 
recently revisited using comparison parabolic arguments in [2] (including the 
case of compactly supported initial data), and [T| (for slowly decaying initial 
data) . Roughly speaking, for compactly supported initial data with convex and 
bounded support, as e — > 0, the solution of (|1.5[) generates a sharp interface at 
the very early stages of the dynamics. Then the interface propagates through 
the spatial domain, according to a free boundary problem with constant speed 
in the normal direction. This speed turns out to be the minimal speed of prop- 
agation of some underlying travelling wave solutions. 

In the delayed case (r > 0) that we consider, we will show that the above 
scenario remains valid under the following assumption on the initial data ip 
arising in (|1.2I) . 

Assumption 1.1. We assume that if : [— r, 0] x M. N — >• [0,1] is a uniformly 
continuous function satisfying the following. 
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(i) There exists w a G BUC 2 (R N , K) sucft i/iai 

ft := {x G R N : w (x) > 0} 
is a nonempty smooth bounded and convex domain, and 

w (x) < <p(0,x), V(0,x) g [-r,0] x M w . (1.6) 

(mJ There exists 5 > smc/i i/iai 

|Vitfo(s).i/ 8 n„(af)| > Vx g T := dft , (1.7) 
wherein Vqq (x) denotes the outward unit normal vector to Oo at x G IV 
(m) There exists vq g BUC(M. N , [0, 1)) smc/i that 

supp u = O , (1.8) 

and 

<p(6,x) < v (x), V(0,x) G [-r,0] x R N . (1.9) 

Remark 1.2. The hypothesis ||i>o||oo < 1 hi (Hi) shall be used in the construction 
of upper barriers in Scction[5] Nevertheless, wlien H^olfoo — 

1, our main result 
remains valid under the additional assumption that / satisfies 

f(K u)<K f(u), Vug [0,1], (1.10) 

for some constant A'o > 1. See Remark 15.21 for details. 

Before stating our main convergence result let us give some notations. Under 
assumption (|1.4j) . we denote by c* > the minimal speed of the underlying 
delayed travelling waves (see Lemma 12.31 for details). In particular, there is 
(U*,c*) G C 2 (]R) x (0, co) such that U* is nonincrcasing and 

f (U*)"{z) + c*(U*)'(z) + f (U*(z + c*t)) - U*{z) = 0, Vz G R, 
\U*(~oo) = 1 and U*(oo) = 0. 

Next, for c > 0, we denote by r c := Ut>o(M x ^t) the smooth solution of the 
free boundary problem (see subsection 14 II for details) 



(P c ) 



V = c on T£ 
r tlt=o = r °> 



with V the normal velocity of in the exterior direction, and To the initial 
interface defined in (|1.7[) . Also, we denote by J7( the region enclosed by the 
hypcrsurfacc V[. 

Here is the main result of the present paper (see subsection 12 . 1 1 for the well- 
posedness of the initial value problem (|l.ip - (|1.2p ). 

Theorem 1.3 (Convergence to a propagating interface). Let the nonlinearity 
f be as in (|1.4j) . Let the initial data ip satisfy Assumvtion l 1 . 1{ For each e > 0, 
let u e : [-er,oo) x R N — > K be the solution of (TO ]) -(TO ]) . Then the following 
convergence results hold. 
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(i) For each c G (0, c*) and each to > 0, we have 

lim sup sup |1 — u £ (t, x)\ =0. 

fjij -For eac/i c > c* and eac/i <o > 0, we have 

lim sup sup u £ (t, x)=0. 

A first step towards Theorcm ll.3l consists in proving that, as e —> 0, the initial 
value problem (|1.1[> — (|1.2p generates a sharp interface after a very small time of 
order O (e| In e [ ) . Then, to analyze the propagation of the interface, we aim at 
constructing suitable sub- and super-solutions. This step is strongly related to 
the existence of travelling wave solutions. While the upper barriers are directly 
constructed by using monostable travelling fronts, the construction of lower 
barriers is much more delicate. This kind of problem has been solved in several 
situations. In [1 1 j . the authors consider a degenerate reaction-diffusion equation, 
and take advantage of the existence of sharp travelling fronts to construct sub- 
solutions. In [2], the standard Fisher-KPP case is considered. The construction 
of lower barriers of propagation is performed by using the existence of non- 
monotone (and also not everywhere positive) travelling waves with speeds c < 
c*. In the non delayed case, the existence of such a connection easily follows 
from a phase plane analysis. In the delayed case we consider, the existence 
of similar waves is far from obvious. The key idea of the present paper is to 
construct sub-solutions of propagation by using travelling waves for a modified 
time delayed reaction-diffusion equation with a bistable dynamics. We hope 
that such a strategy could be used to understand better the classical non delayed 
Fisher-KPP case and also to analyze a larger class of equations. 

The organization of the present paper is as follows. In Section [5J we recall 
known facts on the well-posedness of the initial value problem (|l.ip - (|1.2p . We 
also discuss the links between monostable travelling waves associated with /, 
and bistable ones associated with approximations f n of /. This is necessary 
to develop the key strategy mentioned above. In Section [3l we investigate the 
generation of a sharp interface in the very early stages of the dynamics. This 
is strongly related with the underlying delay differential equation. Section U is 
concerned with the study of the propagation of interface from below. We shall 
construct accurate lower barriers by using a bistable approximation. As a result 
of Sections |31 and HJ we shall prove Theorem 11.31 (i). Section [5] deals with the 
construction of upper barriers to control the propagation from above. This will 
imply Theorem II .31 (ii). 

2 Preliminary 

2.1 Existence and comparison for fll.ip — f |L2|) 

We first state the following comparison principle for monotone delayed reaction- 
diffusion equations. 
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Proposition 2.1 (Comparison principle). Let r > 0, T > and g : R — ► R an 

increasing and continuous function be given. Let (u,v) € C([— t, T] x R^) &e 
too bounded functions satisfying 

d t u, d t v, Vu, Vu, £> 2 u, L> 2 v e L 2 oc ((0,T) x R w ) . 

(9t - A + 1) u(t, x) - 5 - r, a)) < 
(9t - A + 1) v(t, x)-g (v(t - r, ar)) > 0, 

/or almost every (t,x) G (0,T) x M. N , and 

u(6, x) < v(6, x) for all (6, x) G [-r, 0] x R w . 

TTiera u(t,a;) < /or all (t, x) G [— r, T] x R N . 

Proof. Let us consider the map w := u — v G C ([— r, T] x R w ) 
increasing, it follows from (|2.ip and (|2.2I) that w satisfies 

(d t - A + 1) x) < a.e. in (0, min(T, r)) x R N . 

Since w(0, •) < 0, the weak comparison principle [H] Proposition 52.8] implies 
w < in (0, min(T, r)) x M. N . If T > r, one can repeat the argument on 
(r, min(T, 2t)) x R n . This proves the proposition. □ 

We now introduce some notations. Let X := BUC^R^jR) be the Ba- 
nach space of bounded and uniformly continuous functions from M. N to R, en- 
dowed with the usual suprcmum norm. Define also the Banach spaces C := 
C ([— r, 0], X) and Cq := C ([— r, 0], R). For convenience, we identify ^ G C as a 
function from [— r, 0] x R w into R defined by x) = ij){6)(x). For each a < f3, 
we define 

[ aj 0] c := {V GC : a < < p, V(6»,x) G [-r, 0] xR N }, 

and [a,/3]c := Co PI [a,/3]c- Next, for any continuous function w : [— r, oo) x 
R N -> R, we define w t G C, t > 0, by 

u> t : (0, a:) G [-r, 0] x R N ^ w t (6, x) = w{t + 6, x). 

The wcll-poscdncss of the initial value problem (|l.lj) - (|1.2j) can classically be 
investigated via the theory of abstract functional differential equations: since the 
initial data ip G [0, l]c, the initial value problem (jl.ll) - (jl.2|) admits a unique mild 
solution u € : [0, oo) x R N — >• [0, 1], which is actually classical on [et, oo) x R^. 
For more details, we refer the reader to the monograph of Wu and the 
references cited therein. 

2.2 Monostable and bistable delayed travelling waves 

As explained in the introduction, the construction of lower barriers is far from 
obvious when r > 0. A key idea of the present paper is to derive the monostable 
propagation of the interface from below from the bistable case. To perform this 
in Section 31 let us first define a family of bistable approximations by extending 
the monostable nonlincarity / for negative values of u. 



(2.2) 
Since g is 
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Bistable approximations of /. For 77 £ (0,1], we introduce an increasing 
and bounded map f n : R — > R of the class C 2 such that 

/„(«) = /(«) e [0, l] 

/»,(-??) = -V and fv'(-v) < 1 

/^(u) < u Vu £ (-77, 0) U (1, oo ) and /^(u) > u Vu £ (— oo, —77) U (0, 1). 

(2.3) 

Observe that f v has exactly three fixed points —77 < < 1, /?/(— 7?) < 1 and 
f n '(l) = /'(l) < 1. We also require that the family {/^l^wg 1] ^ s ordered in the 
sense that: 

V(r ?J 7 ? ') £ (0,1] 2 , 77 < 77' =* /„,(u) < /„(u) g M. (2.4) 

Travelling waves. We consider the one dimensional reaction-diffusion equation 
with time delay 

(d t - d xx + 1) u(t, x) = f n (u(t — t,x)), t>0, iei. (2.5) 

We denote by u n = u ri (t,x;tp) : [— r, 00) x R N — > [—77, 1] the solution of (|2.5[) 
with the initial condition 

u (6)(x) = u(9, x)=^e [-77, l] c . (2.6) 

Let us notice that the above initial value problem generates a strongly con- 
tinuous and increasing semiflow {Qr){t)} t>0 defined by 

[Q ri (t)1>](9,x) = (u TI ) t {e,x;il>), (M) £ [-r,0]xl s 

and acting [—77, l]c into itself. Also, it follows from (|2.3[) that, for each t > 0, 
Q,,(t)[0, l]c C [0, l]c and that := (Q v (t)) |ro,ii c does not depend upon 77. 
Note that Q»; exhibits a bistable dynamics while Q is of monostable type. 
Let us state some basic facts on travelling waves sustained by (|2.5[) . 

Lemma 2.2 (Bistable Travelling waves). For 77 £ (0, 1] arbitrary, the following 
holds. 

(i) There exists a unique speed c v such that (|2.5[) has a travelling wave solution 
(U n ,c v ) £ C 2 (R) x R whose profile U v is nonincr easing, that is 

(U v "{z) + c n U v '{z) + f v {U v (z + c,r)) - U v (z) = 0, Vz £ R, 
I U v (—oo) = 1 and U v (oo) = —77. 

(nj There exist two constants (n,M) £ (0, oo) 2 sitc/i that 

|1 - L/„(z)| + I-77- | < Me"*, Vz < 0, 

\U v '(z)\ + \U v "(z)\ < Me-f*W, Vz £ R. 

(Hi) There exists some constant 7 > such that, for any ip £ [—77, l]c luii/i 

liminf min ip(9,x) > and limsup max ip(9,x) < 0, (2-8) 
ee[-T,o] 9e[-T,o] 

one can /md if = K(ip) > and £ = £(t/>) £ R suc/i i/iai 

|7i l; (i,x;i/;) - t/^a: - c v t + < A'e" 7 *, V(t, x) £ [0, 00) x R. 
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Proof. Part (i) comes from Schaaf [THl Theorem 3.13] (see also Fang and Zhao 
[5J Theorem 6.4]). The behavior of the profile (ii) can be found in Hupkes and 
Lunel |12[ Proposition 2.2.5]. Finally the global asymptotic stability with phase 
shift of the wave (Hi) is due to Smith and Zhao |181 Theorem 3.3]. □ 

We recall that / satisfies (|1.4[) . As far as monostable travelling waves sus- 
tained by 

(d t - d xx + 1) u(t, x) = f (u(t — t,x)), t>0, (2.9) 

are concerned, we quote the following result from Schaaf [T|)J Theorem 2.7] (see 
also [13]). 

Lemma 2.3 (Monostable travelling waves). There exists c* > such that (|2.9p 
has a travelling wave solution (U c ,c) E C 2 (R) x (0, oo) with < U c < 1, if and 
only if c > c* . In addition, when c > c* the waves are nonincreasing. 

In the sequel we denote by (U* , c*) the monostable wave with minimal speed, 
that is 

(U*)"(z) + c*(U*)'(z) + f (U*(z + c*t)) - U*(z) = 0, Vz E K, 
U*(-oo) = l and U*{oo) = 0. ' 

To conclude this preliminary, we prove the following result on the conver- 
gence of the bistable speeds c v . 

Lemma 2.4 (Convergence of speeds). Let f satisfy (|1.4[) . Let {/jj}jje(o,il satisfy 
(|2.3p and (|2.4p . T/ien the family {c, ; }, fJ(E ^ ji is decreasing and 

c v /• c* , as r]\0. 

Proof. Let 77 E (0, 1] be given. Since < U* < 1 and fr,\[o,i] = /, U*{x - c*t) 
solves (|2.5p . We can select a tp E [—77, l]c such that (|2.8p holds together with 

V>(0,x) < U*(x - c*6>), V(0,x) E [-r,0] x M w . 

The comparison principle yields u^t, x; t/>) < f/*(x — c*t), so that Lemma 12.21 
(Hi) implies 

U v (x - c v t + £) - Ae" 7 * < U*{x - c*t), 

for some constants 7 > 0, A > and £ £ R. Choosing a; = c*t, we get 
C/^((c* - c n )t + £) - A'e -7 * < E/*(0); if c* < c,, then letting t -» 00, we collect 
1 < t/*(0), a contradiction. Hence, we have c r; < c*. 

Now, let us take 77 < 77' in (0, 1]. In view of (|2.4p . the comparison principle 
implies u v >(t,x;ip) fs ^(t, x; VO for any ?/> E [ — 77, Choosing ^ given by 
t/>(#, x) = [/^(x — c^^) and using Lemma 1 2. 2 1 (Hi), we infer that 

U v >(x-Cr,>t + €') - AV 7 '* < U v (x-Cr,t), 

for some given constants 7' > 0, K' > and (el. Choosing h E K such that 
£V(/i) = 0, x = cyi - £' + h, we get -A'e- 7 '* < U n ((c n > - c n )t - £' + ft)); if 
C77' > then letting t — > 00, we collect < —77, a contradiction. Hence, we 
have c„i < c„. 
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As a result, there is c < c* such that c v /* c, as 77 \ 0. To conclude let 
us make the normalization U n (0) = 1/2 for each r\. Classically, by the interior 
elliptic estimates and Sobolev embedding theorem, we may assume that, modulo 
extraction, U n — > U strongly in C\f c (R) and weakly in Wf^ (R), 1 < p < oo. 
Then (J7, c) satisfies (|2.10j) with c* replaced by c. Lemma [2.31 then enforces 
c > c*. The lemma is proved. □ 



3 Lower barriers for small times 

The goal of this section is to prove that, after a very short time as £ — > 0, 
the solution u £ : [-er, oo) x -> [0, 1] of (fTT|) - ([i~2"|) is very close to 1 in O 
(roughly speaking). Precisely, the following holds. 

Proposition 3.1 (Generation of interface from below). Let the initial data ip 
satisfy Assumption ] 1 . 1\ (i) — (ii). Denote by d(0,x) the smooth cut-off signed 
distance function to Tq as defined in subsection ^. 1\ (in particular, d(0,x) < 
if and only if x £ fioj- Then there exist Sq > 0, ao > 0, po > and Eq > such 
that, for all e € (0,£o) and all (9,x) G [— r, 0] x M. N , the following holds. 

If d{0,x) < -8 s\]ne\ then 1 - e p " < u £ (a e\ lne| + er + eO, x) < 1. 

The proof shall be given in the end of the present section. The idea is to 
construct a sub-solution based upon the delay differential equation obtained by 
neglecting diffusion in (11.11) . 

3.1 A delay differential equation 

Let us consider the delay differential equation 

^(t)=f(v t (-r))-v(t), t>0, (31) 
^(0=^(0 6 [0,11c, 

where / satisfies (|1.4[) (recall that Co = C ([— r, 0], R)). Because of the afore- 
mentioned reason, we also need to consider, for ?] e (0, 1], the delay differential 
equation 

dv 

-^(t)=fv(M-r))-v(t), t>0, 

M-) = t K-) e h»7.1]co. 

where was defined in (|2.3p . From standard results for delay differential equa- 
tion with quasi-monotone nonlinearity — see for instance the monograph of 
Smith [T7] — the following holds. 

Lemma 3.2 (Wcll-poscdness). For each <fi e Co, (|3.2j) has a unique global (mild) 
solution Vrf = Vr](-;^)) : [— r, oo) — > R and the semifiow V n (t)(j) — V 71 {t](j)) :~ 
(v n )t(- ; 0) is strongly continuous and monotone increasing on Co- It furthermore 
satisfies the following properties. 

(i) For each t > 0, V v {t)[-i h l] Co C [-T], l] Co . 
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(ii) For each t > 0, Vr,(t)[0, l] Cn C [0, l] Co . The restriction V(t) := V v (t)\ [0Aka 
does not depend upon r\ and, for <fi £ [0, l]e , the map t i-> V(t)(f> = V(t; <j>) 
is the mild solution v t (- ;4>) °f P-1P - 

Dynamics of the DDE. We start with a lemma on the global dynamics of 
(El) on [0,l] Co . 

Lemma 3.3 (Stability of 1). The following holds. 

(i) For (p £ [0, l]c \ {0}, we have lim t _ i>00 V(t)<f> = 1 in Co- 
(ii) There exist S± > 0, M > and A > such that, for all <fi £ Co, 

||W||l~(-t,o)<*i => Hl-^(*)0lU~(-r,o) < Me~ xt , Vi>0. 

Proof. Let us prove (i), that is the global stability of the stationary point v = 1. 
First, we consider the case where there is £ € (0, 1) such that (f)(0) > £, for all 
6 £ [— t, 0]. Since the semiflow associated with (|3.1[) is monotone increasing 
and since V(t)[0,l]c a C [0, l]c , it is enough to consider the solution with the 
constant £ as initial data, that is V(t;() = v t (-;(). Since /(C) > Cj the map 
t <— > v(t;Q is nondecreasing. Hence we get limt_ ! . 00 v(t; £) = 1, which in turn 
implies ||V^(i)C — l||oo = su P- r <e<o \ v (t + ^> — 1| — > 0, as i — > oo. This 
concludes the proof of (i) for this first case. Let us now consider the general 
case. Since (p £ [0, l]c \ {0}, there exist — r < a < < and /3 > such that 

<£(0) >/31 [o , 6] (0), [-T,0]. 

From p. II) . we obtain that, for all t £ (0,r], 

J f (e*v(t;<f,)) > e*/(^)l[r+a,r+6](t) > /(/3)l [T+a , T+ 6](t). 

Integrating this from to r yields v(t; (f>) > f(f3)(b — a). Now, for all t £ (r, 2r], 
p~Tj) implies ^ (e*v(t;0)) > 0. Hence 

«(*; 0) > e T ~MT; 0) > C := e~ T f((3)(b - a) > 0, Vt G [r, 2r], 

and we are back to the first case. This completes the proof of (i). 

The proof of (ii) is a direct consequence of the exponential local stability of 
v. Indeed, at this point the characteristic equation associated to (|3.1[) reads as 

A (A) := A + l-/'(l)e- Ar = 0. 

Since f'(l) < 1, all roots have strictly negative real parts and the result follows 
(see for instance [3T], [TU] and the references therein). □ 

Next, we shall prove the following important result. 

Proposition 3.4 (Convergence to 1). Let <j> > in Cq \ {0} be given. There 
exists A > such that, for all a > there exists Eq = £o(oe) > such that, for 
all e £ (0, £o), 

1_ £ «V2 < V (a\\ne\+t;e\]ne\<t>) (0) < 1, V(0,i) £ [-r,0] x [0,oo). 
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Proof. Let > in C \ {0} be given. Recalling that /'(0) > 1, let 6 e (0, 1) 
and p > 1 be such that 

f(u) > pu, Vw € [0,(5]. (3.3) 

Applying Lemma 13.31 with 5 as initial data, we have the existence of constants 
M > and A > such that 

< 1 - V(t; 5)(6) < Me- Xt , V((9, t) e [-r, 0] x [0, oo). (3.4) 

Let a > be given. Consider £o > small enough so that e\ lne|^> £ [0, 8]c for 
all e £ (0, £ ). Since > is in C \ {0}, there exist -r < a < 6 < and /3 > 
such that 

e\lne\4>(e) >s\ln £ \pi [aM (e), V6 € [-r,0]. 

Arguing as in the proof of Lemma 13.31 and using (|3.3[) . we discover that there is 
C > such that, for e > small enough, 

v e (t) :=v(t;e\1ns\<f>)>Cs\]ne\, Vie[r,2r]. (3.5) 

Next, observe that, for all < t < t, 

j t (e*v e (t)) = e t f{e\ \ne\4>(t - t)) < e T e\ lneH^IUH/'Hoc =: Ce\ lne|. 

Integrating this from to r, we have v e (r) < e~ T ((j)(0) + Cr)e\\ne\ < 5, for 
£ > small enough. Therefore we can define 

t £ := sup {t > 2r : v e (s - r) < 5, Vs E [2r, t]} . 

It then follows from the DDE (gH]) and ([23]) that 

u £ '(t) > pu £ (f -r) - u e (i), Vt e [2r,t e ]. (3.6) 

Since p > 1, there is a > such that a + 1 = pe~ aT . Then the map h : t ^ 
Ae\ lne|e a *, A := C/e 2or satisfies 

h'(i) = ph(t-T)-h(t), £ [2T,i e ], and /i(t) < C£|lne|, te[r,2r]. (3.7) 

It follows from (O, (j33J) and (jHUl) that 

u e (i) > Ae|lne|e Qt , Vt £ [2r,t e ]. 

In view of v £ (t e — t) = 6, we have 

1 5 

t £ <r + -In— r. (3.8) 

~ a Ae\\ne\ 

Now since the map t t-y v e (t) is increasing, we deduce from v e (t £ — t) = S that 

v s (t e +t + 6)>8, V(6», t) e [-t, 0] x [0, oo). 
In view of (|3.8[) , we have t e < a \ In e \ for e > small enough so that 
v £ {a\\ns\+t + 9) > S, V(6»,t) £ [-r,0] x [0,oo). 
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Since the semiflow associated with (|3.1|) is monotone increasing on Co, we thus 
have 

< l-v e (a\1ne\+t + e) < 1 - V(a\ \ne\ + t; 5){d), 
which combined with Q3.4p yields, for e > small enough, 

< 1 -v s (a\\ns\ +t + 6) < Me~ A(a|llie|+t) < Me aX < e aX/2 . 

This completes the proof of Proposition 13.41 □ 

Derivatives of the semiflow. Let us now provide some estimates on the 
derivatives of the semiflow with respect to the state variable. Our first result 
is a consequence of the well-known differentiability result of semiflows generated 
by delay differential equations (see for instance [TIT], see also [2TJ for results on 
abstract scmilincar problems with Hille-Yosida non-densely defined operator). 

Lemma 3.5 (Derivatives). For each t > 0, the map <f> £ Co ^ V v (t;<fi) £ Co 
provided by Lemma VS '.2\ is of the class C 2 . For each <fio £ Co and each £ Co, the 
map t £ [0, oo) i — ^ d^V^t; 4>q)-4> £ Co is the mild solution of the non- autonomous 
equation 

dv 

-(t)=L(t, (f> o)v u t>0, (3g) 

v{e) = m, 0e[-T,o], 

wherein, for each t > 0, L{t, <po) : Co — > M is defined by 

L(t, ^ := (V v (t; o )(-r)) 0(-r) - 0(0). (3.10) 

Moreover, for each <fio £ Co and each <f> £ Co, the map t ^ d 2 ^V^t; 4>q) ■ (0, <f>) 
is the solution of 

dv 

— (t) =L(t,0oH + G(t;0 o ;0), t>0, (3U) 
v(9) = 0, 6 £ [-r,0], 

wherein the map 1 1— > G(t; </>o; 0) *s defined by 

G(t; O ; 0) := /,/' (V„(t; fo)(-r)) [0*V„(f; 0o) • 0(-r)] 2 . (3.12) 

Here is an estimate on the first derivative. 

Lemma 3.6 (First derivative). There exist constants M + > 1 and j + > such 
that, for all 4>q £ Co , 

e -r e -(t+e) < d^V v {t-fa) ■ l{6) < M + e^ + ^ t+e \ V(0,t) £ [-r,0] x [0,oo). 

Proof. Let 0o £ Co be given. First, the semiflow V^(i) being monotone increas- 
ing on Co, observe that 

d^,V v (t; O ) ■ 1(6>) > 0, V(0, t) £ [-r, 0] x [0, oo). (3.13) 

Hence, in view of (|3.9[) and (|3.10j) . the function w(t) := d^,V v (t; O )T(O) satisfies 

w'(t) > -w(t), Vi > 0, 
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so that w(t) > e *, for alH > 0, which in turn implies 

diV v (t;<M-l(0) >e- (t+9) , 

for all (9,t) £ [-r, 0] x [0, oo) such that t + 9 > 0. For the remaining (0,t) £ 
[-r,0] x [0, oo) such that t + 9 < 0, we have fyV^i; (f> a ) ■ 1(9) = 1 > e -( T +*+ e ). 
This completes the proof of the left-hand side of the estimate of the lemma. 
Next, choosing a constant N > 1 such that 

< f v '(u) < N, Vu £ M, (3.14) 

we infer from (|3.9|) and p.lOp that 

w'(t) <Nw(t-r)-w(t), t > 0, and w(6) = 1, 0e[-r,O]. (3.15) 

Observe that the map h :t i-> e^ JV_1 ^ T e^ JV ~ 1 ^ satisfies 

h'(t)>Nh(t-T)-h(t), t > 0, and > 1, fle[-r,0]. (3.16) 

It follows from (j3~T5|) and (|3~T6)) that w(t) < eC^- 1 )"^^- 1 )*, for all t > . 
Arguing as above we get the right-hand side of the estimate of the lemma. □ 

We pursue with the following estimate on the second derivative. 

Lemma 3.7 (Second derivative). There exist constants K > and /i > such 
that, for all (fto £ Cq, 

\d^V v (t;cf> ) ■ (1,1)(9)\ < Ke^ t+e \ V(0,i) £ [-r,0] x [0,oo). 

Proof. In view of (|3.12[) and Lemma 13.61 there exists a constant A > such 
that, for all cf>o £ Cq, 

|G(t;0 o ;l)| < Ae 27+(t - r) , V< > 0. 
Hence, the function w(t) := d^V v (t; cj>o) ■ (1, 1)(0) satisfies 

w'{t) < Nw(t - r) - w(t) + Ae 27+(t " T) , i > 0, and u;(0) = 0, 9 £ [-r, 0]. 

(3.17) 

We look for a super-solution of (|3.17[) in the form 1 1-> i^e' 14 , for some constants 
K > and /i > to be determined. This leads us to 

/ t >A>e-' iT -l + 4-e- 27+r+ ( 27+ -' i )*, Vt>0, (3.18) 
K 

which can be achieved by choosing fi > 2j + and K > both large enough. 
Arguing as in the proof of Lemma 13.61 we end up with constants K > and 
yu > such that, for all </> £ Co, all e [— r, 0], all f > 0, 

a^K ) (t;0 o )-(i,i)(0)<^ t + e ). 

Next, select C > such that f v "(u) > —C, for all u G M. Then we get 
w'(t) > —Cw(t — t) — w(t) — Ae 2l+ ^~ T \ for which we can construct a sub- 
solution t h-> —Ke^ as above. This completes the proof of the lemma. □ 
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As a direct consequence of Lemma 13.61 and Lemma 13.71 we obtain the fol- 
lowing estimate. 

Proposition 3.8 (Estimate on derivatives). There exist constants K > and 
7 > such that, for all 4>o G Co, 

\d^V v (t; 4> ) ■ (1, 1)(0)| < Ke^V^t; O ) ■ 1(0), 

/or (0,i) e [-t,0] x [0,oo). 

3.2 Construction of lower barriers for small times 

We now provide an accurate lower estimate, for small times, of u E : [— er, oo) x 
R N -> [0, 1] the solution of (p~Tj) - (fT2|) . 



Proposition 3.9 (Sub-solutions). Lei the initial data ip satisfy Assumvtion \l.l\ 
(i). Then there exist K > 0, a > and Eq > smc/i i/iai, /or aZ/ e £ (0,£o), 

max |o ; v r) f-; lOoO^) — ei^r — iff J } < u £ {t, x), 

for all (t,x) £ [— er, ae\ lne|] x IR^. flere, u,, = v v (-;<fi) : [— r, oo) — » R denotes 
the solution of (|3.2p arising in Lemma \3.S\ and the function wq is as in (jl.6p . 

Proof. Let us consider the differential operator 

£^ [it] (i, x) := d t u(t, x) — eAu(t, x) — - [/ r; (w(i — er, x)) — it(t, x)] . 

Since f v = f on [0, 1], we have [u s ] (i, x) = 0. We look for a sub-solution, at 
least for small times, u : [—er, oo) x M. — > K in the form 

u(t, x) := v v ( — ; u)q(x) — eKt — Kt 



Straightforward computations yield, for each t > and each x S 



C e v [u](t,x) = -V £ (t,x) 



W £ (t t) 

K + eAw (x) + e l ' ! \Vw (x)\ 2 
V e {t,x) 



1 

+ - 



dv„ \ /i , , „ \ . ( ft 



; iy (x) - et - Kt - f v [v v [ r;w (x) - Kt 



where 



V e {t, x) := d^V v (pw (x) - eKt - Kt 
W 6 (t,x) := d^V v (^;w (x) - eKt - Is 



Since the semiflow arising in Lemma 13.21 is monotone increasing in Co and since 
fq is increasing, we have 

dv v \ ft , s „.\ , f ft 



dt 

0, 



w (x) - eKt - Kt - f v [v n \ t; w (x) - Kt 



< ( + v v ) ( ~; iu (x) - eKt - Kt] - f v (v v [ - - r; iy (a:) - - 
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since v v solves (|3.2p . Hence, using Proposition 13.81 we get, for all e £ (0,1), 

(>0,ieK w , 



C e v \u](t,x)<-V E (t,x) 



A--e||A U ;o||oo-e||V U ;o||L^ 



Looking at small times, the above implies, for all e £ (0, 1), t £ (0, 7 1 e|ln£|), 



x £ 1& N , 



C £ v [u](t,x) < -V e (t,x) K - e\\Aw \\oo - WVwoW^K 



< 0, 



if K > is sufficiently large. Next, concerning initial data, we have, for all 
0e[-er,O], 



u(0, x) = wq{x) — eKt — K6 < wq(x) < ip I — , x j = u e (9, x), 

where we have used p.6|) and (|1.2p . The comparison principle in Proposition 
12.11 thus implies that 

u[t,x) < u e (t,x), V(t,x) £ [— er, 7 _1 e| lne|] x R N . 

Recalling that u e > 0, this completes the proof of Proposition 13. 91 □ 

Proof of Proposition 13. ll Fix K > and a > as in Proposition ^. 91 Define 
Q'o := a/2. For <j> := ao £ Co \ {0}, let us select A > as in Proposition 13. 41 and 
define po := aoX/2. Also, it follows from Assumption 11.11 (ii) that there exists 
Sq > such that, for e > small enough, 

d(0,x) < -S e\lne\ =^ w Q {x) > 4a e|lne|. (3.19) 

Now, for any — r < 8 < 0, define s := aoe|lne| + er + eO and take x such 
that d(Q,x) < — 6oe\lne\. Since, for e > small enough, < s < ae|lne| and 
wo(x) — eKt — Ks > aoe \ lne|, we deduce from Proposition l3.9l and Proposition 
[531 that 

u 6 (s, x) > v v (a \ lne| + t + 9; a s\ lne|) > 1 — e p " , 
which concludes the proof. □ 

4 Lower barriers via bistable approximation 

As explained before, our analysis of the propagation of interface from below is 
performed by approximating the monostable function / in a bistable manner (see 
subsection I2.2[) . We start with some preliminaries on smooth signed distance 
functions associated with a family of free boundary problems. 

4.1 Smooth cut-off signed distance functions 

For c > 0, we denote by r c := Ut>o(M x ^t) * ne smooth solution of the free 
boundary problem 



(P c ) 



V = c on 
r * lt=o = r °' 
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where V denotes the normal velocity of in the exterior direction. Note that 
since the region enclosed by To, namely £Iq, is convex, these solutions do exist 
for all t > 0. Also we can naturally, i.e. in a reversible manner, extend these 
solutions for small negative times by letting r evolve with speed — c. Hence, 
with a slight abuse of notation, we consider for all t > —st, with e > 
small enough. For each t > —st, we denote by fij the region enclosed by the 
hypcrsurfacc T^. 

Let d be the signed distance function to T c defined by 



-dist(x,r^) forarenj 
dist(>,r£) for x eR N \ n r t , 



where dist(a;, r£) is the distance from x to the hypersurface T^. We remark that 
d = on r c and that |Vd| = 1 in a neighborhood of T c . 

We now introduce the "cut-off signed distance function" d, which is defined 
as follows. Let T > be given. First, choose do > small enough so that d is 
smooth in the tubular neighborhood of F c 

{(t,x) G [-er,T] x R N : \d(t,x)\ < 3d }. 

Next let £(s) be a smooth increasing function on R such that 

{s if \s\ < do 

-2d if s < -2d 
2d if s > 2d . 

We then define the cut-off signed distance function d by 

d(t,x):=c(d(t,x)y (4.2) 

Note that 

if |d(t,a:)| < «f then |Vef(i,a;)| = 1, (4.3) 
and that the equation of motion (P c ) yields 

if \d(t,x)\<d then d t d(t, x) + c = 0. (4.4) 

Then the mean value theorem provides a constant N > such that 

\d t d(t,x) + c\ < N\d(t,x)\ for all (t,x) G [st,T] x R N . (4.5) 
Moreover, there exists a constant C > such that 

\Vd(t,x)\ + \Ad(t,x)\ <C for all G [-er,T] x R w . (4.6) 

4.2 Construction of lower barriers 

Let us recall that {/rj}r7£(o,i] denotes a family of bistable approximations of / 
such that (|2.3p and (|2.4p hold. Also, for r\ G (0, 1], (L^, c v ) denotes the travelling 
wave solution (with time delay) associated with this bistable (see Lemma l2~2"]) . 
namely 

f C/ 1) "(z) + Ci; C/ 1) '(z) + / t? (C/ 1) (z + c 1) r))-C/ I) (z) = 0, Vz G R, 

1 E/„(-oo) - 1, U n {0) = 0, l/,(c») = -v- [ ' ' 
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In the spirit of the sub-solutions constructed in [5] for bistable systems, we 
look for sub-solutions u~ in the form 

u-(t, X ) := U V ^A±phMi y q{t) , (4 . 8) 

where 

p{t) := -e- fjt ' £ + e Lt + K, (4.9) 

q(t) := a^(3e-^/ £ +eLe Lt y (4.10) 

Here, <r, /3, L and K are positive constants to be determined, and d n (t, x) denotes 
the cut-off signed distance function to the interface starting from To and evolving 
with speed c v , that is the solution of (P Cr >). As seen in the previous subsection, 
this allows to define u~ for all t > — er, x 6 



Proposition 4.1 (Sub-solutions). One can find positive constants (3, a and L 
such that, for all K > 1, the function u~ satisfies, for e > small enough, 

e£^[u~](t,x) = ed t u~(t,x) - e 2 Au~ (i, x) - f v (u~(t - er,x)) +u~(t,x) < 0, 

for all t > 0, x £ R N . 

Proof. For ease of notation, we drop most of the subscripts rj. Also we define 

z ._ d(t,x)+e\\i\e\p(t) . 

£ 

We start by evaluating eC^[u~](t,x). We compute 

ed t u~{t,x) = (d t d(t,x) + e\ lne\p'(t))U'{z) - eq'{t) 
e 2 Au-(t,x) = \Vd\ 2 (t,x)U"(z) +eAd(t,x)U'(z). 

Next, observe that the previous subsection enables to write 

d(t — er, x) = d(t, x) + ecr + £0 E (i, x), 

where the correction 9 e vanishes close to the interface and is 0(1): 

9 e (t,x)=0 if \d(t,x)\ <d , ||e e ||ioo<A (4.12) 

for some constant A > 0. Hence, since p(t) increases and U (z) decreases, we 
have 

, . TT { d(t, x) + e\ lnelpft — er) _ , ,\ , . 

u-(t-er,x) = U[-^—^ ■ — -+CT + e E (t,x) ) - q(t - er) 



TT { d(t, x) + e\ \ne\p(t) ^ . ,\ , . 

> rj f _!i2_J — J +CT + 0£ ( t; x ) \ -q(t- er). 

Since / is increasing we get 

f (u-(t-ST,x)) > f(u(z + cr + Q e (t,x))-q(t-erj) 

= f(u{z + CT + Q e (t,x))) -q(t-sr)f'(d), 
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for some U(z + ct + Q e (t, x)) — q(t — et) < 9 < U(z + ct). Hence, we have 

f(u-(t-er,x)) > f(U(z + cT))-q(t-er)f'(9) 

+Q £ {t, x)(f o U)' (z + ct + uQ e (t, x)) , 

for some < co < 1. Combining the above estimates with U"(z) + cU'(z) + 
f(U(z + cr)) - U{z) = 0, we obtain eC^[vr](t, x) < E l + E 2 + E 3 where 

E 1 := e\\iie\p'{t)U l {z) + q{t~eT)f'{e)-q{t)-eq l {t) 

E 2 := {d t d{t 1 x)+c~eM{t 1 x))U\z)+(\-\Vd{t 1 x)\ 2 )U"{z) 

E 3 := -Q e (t,x)(foU)'(z + cr + ujQ £ (t,x)). 

Let us now analyze further the term E\. By using the expressions (|4.9D . 
(|4.10p for p and q we obtain 

E x = pe-W e (\lne\U'(z) + *(eP T f'(9)-l + (3)) 

+eLe Lt (| \ne\U'{z) + a( e - eLr f \6) - 1 - eL)) 
=: Pe-^^h+eLe^h. 

Since f'{—rf) < 1 and /'(l) < 1, we can fix small a > and j3 > such that 

e 0T f'{u) -l + p<-p, Vue [-77 - a, -r? + a] U [1 - a, 1 + a]. 

In view of £7(— oo) = 1, U (oo) = — 77 and inequality U(z + ct + Q s (t, x)) — q(t — 
et) < 9 < U(z + ct), there exists a large zq such that 9 £ [—77 — a, —77 + a] U [1 — 
a, l + a] as soon as |z| > zq (by choosing a small enough to control the — q(t — et) 
term) and the above inequality applies for s = 9. It follows from U'(z) < that 
I\ < —a[3 in the region {\z\ > zq}. In the compact region {|z| < zq}, we have 
U'{z) < —6 for some b > so that 1\ < —b\ lne| + C so that ii < — aft also 
holds true. The same argument yields I2 < — cr/3. Hence 

Si < -cr/3 2 e- /3t/e - eaPLe Lt < -sa/3L. 

We now conclude the proof of eCt [u~)(t, x) < 0. Assume first that (t,x) 
lies in the tubular neighborhood {\d(t, x)\ < do} of T t . In view of (|4. 3[) and 
(|4.4|) . the term i? 2 reduces to -eAd(t,x)U'(z). In view of (|4. 12f) . the term E 3 
vanishes. As a result, 

e£*[u-](t,x) < -eapL + £\\Ad\\ L eo\\U'\\ L oo (R) < 0, 

if L > is large enough. Next, if (t,x) is such that \d(t, x)\ > do then we 
shall use the exponential decay of the derivatives of U — see Lemma l!?T21 (ii) — 
to control E 2 and E3. Indeed in this region, the argument z defined in (|4.11[) 
satisfies \z\ > do /(2s). Hence, combining the exponential decay of U' and U" 
with (|4~5|) and (|4.6p . we get a bound \E 2 \ < C 2 e~ C2 for some C 2 > 0. Also, 
it follows from P~T2|t that 

|z + ct + u>Q e (t, x)\ > ^ — CT — £t) ^4 > J^, 



17 



which in turn provides a bound I-E3I < C^e C3 *= , for some C3 > 0. As a result 
we collect, for a constant C > 0, 

e£* [«-](*, a;) < -ecr/?L + Ce~ (7 ^ < 0, 

if £ > is small enough. This completes the proof of the lemma. □ 

In order to apply the comparison principle, we need the following estimate. 

Lemma 4.2 (Ordering initial data). One can find K > 1 such that, for e > 
small enough, 

u~(t,x) < u £ (t + a e\ lne| + et,x), for all - et < t < 0, x £ M JV , 

where aos\ lne| denotes the "generation of interface from below time" appearing 
in Proposition \ 3. 1[ 

Proof. For ease of notation, wc drop most of the subscripts n. If (t, x) is such 
that d(t, x) > —e\ \ne\p(t), then the decrease of the wave U yields u~(t, x) < 0, 
and there is nothing to prove. Now let us take (t,x), with — et < t < and 
d(t,x) < — e\hie\p(t). From the generation of interface from below analysis we 
know that (see Proposition 13. ip 

d(0,x) < -6 e\\ns\ ==^ 1 - e Po < it e (a e|lne| +er + t,x) for - et < t < 0. 

(4.13) 

Writing d(0, x) = d(t, x) + 0(t) and using the expression for p in (|4.9[) . we get, 
for — et < t < 0, 

d(0,x) < -e\ In e\p{t) + Cer 

< -e\hie\(-eP T + e~ eLr + K) + Cer 

< -S e\lne\, 

for £ > small enough, if K is chosen sufficiently large. In view of (|4.13|) it 
suffices to show that u~ (t, x) < 1 — e Po , which follows from the vertical shift q. 
Indeed, the expression for q in (|4.10[) shows that q(t) > a (3 for — et < t < 0, so 
that u~(t, x) < 1 — a/3 < 1 — e p " . The lemma is proved. □ 

Proof of Theorem 11.31 (i). From Proposition 14.11 Lemma T4.2I and the com- 
parison principle, we infer that 

u~(t - a e\ lne| - et, x) < u £ (t, x) for all f > a e\ lne| + et, x G R N . (4.14) 

Let us recall that u~ is defined in (|4.8p and that U n (—oo) = 1. Hence, the 
convergence to 1 in fi£ , as expressed in Theorem[0](i), is a direct consequence 
of both Lemma [^41 and the lower estimate (|4.14p . □ 

5 Global in time upper barriers 

The aim of this section is to construct a super-solution in order to control 
the propagation of the solution from above. Let (U*,c*) be the monostable 
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travelling wave with the minimal speed c* > (see Lemma 12. 3p , namely 

' (U*)"(z) + c*(U*)'(z) + f (U*(z + c*t)) - U*(z) = 0, Vz G R, 
(U*)'{z) < 0, Vz G R, 
C/*(-oo) = 1 and [/*(oo) = 0. 

Then we shall prove the upper estimate on u E : [—er, oo) x R w — > [0, 1] the 
solution of (fL~T1) - ([L~2")) . 

Proposition 5.1 (Super-solutions). Let the initial data ip satisfy Assumption 
Denote by d(0, x) the smooth cut-off signed distance function to Tq as 
defined in subsection \4- 1\ (in particular, d(0,x) < if and only if x G Qq)- Then 
there exists h G R such that, for all e > small enough, 

u £ (t, x) < U* ^(M~ c * f + ^ t V(tj 3.) e [_ £T5 oo) x R w. 

Proof. Since the function vq appearing in Assumption !!. II fm) satisfies ||wo||oo < 
1, we can choose ligE such that ||vo||oo < U*(c*t + h). Up to changing U* by 
[/*(• + h), we can assume h = so that 

||«o||oo < Z7*(c*r). (5.1) 

Let Xq G c%lo = To be given and denote by no the outward unit normal vector 
to Tq at xq. Then consider the map u + : [—er, oo) x M. N — > R defined by 



o ... (a: - x ).n - c*i 

Setting z = , we compute 

£ 

C £ \u + ](t,x) := d t u + (t,x) - eAu + (t,x) - -f (u + (t - er,x)) + -u + (t,x) 

e e 

= -- (f/*)' (z) - - (U*)" (z) -f (U*(z + c*r)) + V(z) 
see e 

-0, 

for all t > 0, x € M. N . Let us now prove that 

«■(*,*) = v (f *) < tr = 

for all (0, x) G [— er, 0] x M. N . In view od Assumption ll.il (Hi) and the decrease 
of U* , it is sufficient to check that 

ifr(s)<ff*( (g ~* o) '"° +cv), VxgR n (5.2) 

When (x — xo).rio < 0, the above inequality follows from (|5.1[) . When (x — 
Xo).Ho > 0, (|1.8[) and the convexity of fio implies vq(x) = and (|5.2|) is clear. 
Hence, it follows from the comparison principle that 

u E (t,x) < U* ( ^ x - x v)-^~ cH ^ ; y M E [_ £T500 ) x M ^ 
for each xq G SJIq. This completes the proof of the proposition. □ 
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Remark 5.2. If ||vo||oo = 1 then, under assumption f| 1 . 10[) of Remark 11.21 we 
have C € [Kou + ](t, x) > 0. Also, normalizing the travelling wave U* by 1 = 
K U*(c*t) and arguing as above, we see that u e (8,x) < K u + (8,x), for all 
(0,x) £ [— £T, 0] x R N . Hence, the comparison principle yields 

u £ {t,x) < K U* ^ x - x o)-™o-c*t ^ ^ y ^ g [ _ £T)0o) x jjW 
for each x G <9£lo- 

Proof of Theorem 11.31 (ii). The convergence to outside Q£ , as expressed 
in Theorem ll.3l (ii), is a direct consequence of the control from above provided 
by Proposition l5.ll □ 
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